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Abstract. We study the transverse-momentum (qr) and rapidity (j/) distributions of the Higgs bo- 
son in perturbative QCD, including the most advanced theoretical information presently available: 
fixed-order perturbation theory at Next-to- Leading Order (NLO) in the large-qr region {qt ~ Mh, 
being Mh the Higgs mass), and soft-gluon resummation at the Next-to-Next-to-Leading Logarith- 
mic accuracy (NNLL) in the small-gr region (qr <^ Mh)- We present numerical results for the 
doubly-differential {qr and y) cross section for the production of a Standard Model Higgs boson 
at the Large Hadron Collider (LHC). 

PACS. 12.38.Bx Perturbative calculations - 12.38.Cy Summation of perturbation theory - 
14.80.Bn Standard-model Higgs bosons 



Introduction 

The gluon fusion process gg ^ H, which proceeds 
through a heavy-quark loop, is the main production 
mechanism for the Standard Model Higgs boson [1] at 
the Large Hadron Collider (LHC) over the full mass 
range 100 GeV < Mh < 1 TeV [2]. As a consequence, 
in the last decade an enormous theoretical effort has 
been devoted to the computation of higher-order per- 
turbative corrections both for this signal and for its 
main backgrounds, in order to achieve the highest pos- 
sible theoretical accuracy. 

A particularly important observable is the doubly- 
differential transverse-momentum (qx) and rapidity (y) 
distribution: a precise knowledge of the Higgs qx- and 
y-spectrum is very important to improve the statis- 
tical significance at hadron colliders by applying, for 
instance, suitable cuts on the jets accompanying the 
Higgs decay products [2l[3]. 

At 0{a'g), the Higgs is produced with vanishing 
transverse-momentum. In order to have non- vanishing 
values of (7t, a recoiling jet is required and thus the 
Leading Order (LO) transverse-momentum distribu- 
tion starts at 0{ag) The Next-to-Leading Order 
(NLO) QCD corrections have been calculated in the 
infinite top mass hmit (Mt oo) [S l IH l lT l llO p il). i.e. 
by using an effective lagrangian directly coupling the 
Higgs to gluons. This approximation has proved to be 
sufficiently accurate provided that Mh < '2Mt and 
qr < Mt 019]. 

It has long been known that, in the small-gr re- 
gion (qx ^ Mh), the presence of large logarithmic 
terms of the form aslog{M^/q^) spoils the conver- 
gence of the perturbative series. These terms originate 



from the emission of soft and coUinear radiation from 
the incoming partons. Since the bulk of the events is 
expected in the small-gr region, an all order summa- 
tion of the logarithmic enhancements is mandatory in 
order to obtain reliable results. 

The technique to perform soft-gluon resummation 
at small transverse-momentum in perturbative QCD 
is well-known [Tl[Il|Il[Tll|Tll[T71|TSl|in] and has been 
applied to the Higgs case up to Next-to-Next-to-Leading 
Logarithmic level (NNLL) [5^. The fixed-order and 
resummed results have eventually to be matched in 
order to prevent possible double-counting of the loga- 
rithmic terms in the intermediate-gT region and thus 
to obtain a uniform theoretical accuracy over the en- 
tire (7T-range. The mat-ching is achieved by taking the 
sum of the two contributions and then subtracting the 
truncation of the resummed term to the same pertur- 
bative order of the fixed-order result. 

In 21, 22J we provided the details of the transverse- 
momentum resummation formalism that we developed 
for the hadroproduction of a general colourless final 
state, and we performed a detailed phenomenological 
study in the case of Higgs boson production at the 
LHC. We included the NNLL resummed result and the 
purely perturbative calculation at NLO, thus reaching 
a uniform theoretical accuracy of 0{ag) over the entire 
gr-range. The formalism has been implemented in the 
publicly available numerical code HqT [|23, . Lately [24] 
we extended the resummation formalism to include ra- 
pidity dependence, thus providing NNLL+NLO accu- 
racy for the fully- differential cross section in qx and y 
at the LHC. 
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The inclusion of rapidity does not change the main 
features of our formahsm: 

— the resummation is performed at the level of the 
partonic cross section with factorization of the par- 
ton distribution functions as in the customary fixed- 
order calculations; 

— the formalism can be applied to any hard-scattering 
process producing a colourless final state accompa- 
nied by an arbitrary and undetected final state; 

— the singular terms are exponentiated in a universal 
(i.e. process-independent) form factor; 

— a constraint of perturbative unitarity imposed on 
the resummed contribution allows both to decrease 
the uncertainty in the matching procedure at inter- 
mediate qx and to recover the total cross section 
result upon integration over qx- 

In the following we will show numerical results for 
the Higgs differential distributions at the LHC. For de- 
tails about the formalism and for further phenomeno- 
logical discussions, we refer the reader to our previous 
papers [2Tl[2l[24]. 



Numerical results 

We present numerical results for the doubly-differential 
{qx and y) cross section for the production of a Stan- 
dard Model Higgs boson with mass Afif=125 GeV at 
the Large Hadron Collider. We used the MRST2004 
NNLO (NLO) set of parton distribution functions [25] 
with as evaluated at 3 loops (2 loops) for predictions 
at NNLL-t-NLO (NLL-I-LO) accuracy. We fixed the 
renormalization and factorization scales both equal to 
the Higgs mass jifj = jip = Mh and made them vary 
between Mh/'2, and 2Mh to examine the scale depen- 
dence of our results. As a cross-check of our calcula- 
tion, we have verified that we reobtain both the nu- 
merical results of Ref. j22j upon integration over y and 
the NNLO total cross section at fixed y [TOlfTT] upon 
integration over qx- 

In Fig.[l]the scale dependence of the NLL-fLO and 
NNLL-I-NLO gr-spectrum with integrated rapidity is 
shown: the reduced thickness of the NNLL-f NLO band 
with respect to the NLL-I-LO one and the overlapping 
of the two bands in the region qx < 100 GeV indicate 
a very good convergence of the resummed result. In 
the upper-right corner we show the K factor defined 

by 

dcTNNLL+NLoif-F, M-R) 



1,0 



Kiqx) 



dcTNLL+LoitJ'F = IJ'R = Mh) ' 



(1) 



i.e., the NNLL-f NLO band normalized to the central 
value of the NLL-I-LO one. We note that a simple 
rescaling of the NLL-I-LO result is not allowed since the 
K factor turns out to be gT-dependent: the cross sec- 
tion is enhanced in the large-qr region, where higher- 
order contributions are not negligible [lllUlTj , and sup- 
pressed in the small-gr region, where the non-per- 
turbative regime sets in. 

In Fig. [2] we plot the qx dependence of the cross 
section at ?;=0, showing both the purely perturbative 
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Fig. 1. Comparison of the NLL-fLO and NNLL-f NLO 
bands. The NNLL-fNLO result normalized to the central 
value of the NLL-fLO result (see Eq. [1]) is shown in the 
inset plot (from [22j). 



NLO result and the resummed NNLL-fNLO result. 
The NLO cross section diverges to — oo as ^ due 
to the large logarithmic terms coming from soft-gluon 
radiation, and shows an unphysical peak. In contrast, 
the NNLL-fNLO result is perfectly regular at small- 
qx, vanishing for qx=^ and converging to the NLO 
result for higher qx values {qx ~ Mh)- The resumma- 
tion effects are clearly visible when looking at the in- 
side plot, where the ratio of the matched NNLL-fNLO 
result to the NLO fixed-order result is shown: 



K(qx,y) 



da 



NNLL+NLO 



/{dqxdy) 



daNLo/{dqxdy) 



(2) 



Resummation is not only relevant at small qx but 
also in the intermediate region {qx <80 GeV), where 
there is a ~20% enhancement with respect to fixed- 
order. The small difference between the y=0 curve 
(solid line) and the integrated rapidity result (dashed 
line) evidentiates the poor rapidity dependence of the 
resummed result. 

In Figure [3] the rapidity dependence of the cross 
section at qx=l5 GeV is shown both at NLO (dashes) 
and NNLL-fNLO (solid line) accuracy. The resummed 
result reduces the cross section in the central rapidity 
region, where most of the events are expected (~25% 
suppression with respect to fixed-order). The K factor 
of Eq. [21 shown in the inset plot, is roughly constant 
in the central rapidity region and starts to be rapidity- 
dependent in the forward (and backward) region where 
the cross section is rather small. This behaviour also 
explains the coincidence of the two curves in the inset 
plot of Fig. [21 

A more detailed investigation of the rapidity de- 
pendence of the cross section is obtained by studying 
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Fig. 2. The Higgs gT-spectrum at the LHC for y—0. The 
K factor defined in Eq. [2] is shown in the inset plot (from 

HU). 
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Fig. 4. The Higgs rescaled gr-spectrum at the LHC as 
defined in Eq. [3] The dotted lines show the values for the 
ratio Ry defined in Eq. |4](from [24]). 
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Fig. 3. The Higgs y-spectrum at the LHC for gT=15 CeV. 
The K factor defined in Eq. [2] is shown in the inset plot 
(from pi]). 



the quantity 



R{qT\v) 



d'^a/{dqT dy) 
da I dqx 



and its QT-integrated version 

da/dy 

Uy = . 



(3) 



(4) 



In Fig. 2] we plot these two quantities, as a function 
of qx, for the two different values y=0 and y—2. The 
NNLL+NLO and NLO results are nearly equal at fixed 
rapidity, reflecting the similar behaviour of the inset 
plot in Fig. [2l The overall decrease of the differential 



cross section when going from y=0 to y—2 amounts to 
nearly 40%, as expected since the total cross section 
rapidly decreases with increasing rapidity. As for the 
qT dependence, the results show a slightly increasing 
(decreasing) slope for y=0 {y=2) and it is quite evident 
that the cross section varies more in absolute value 
than in qx shape. 



Summary 

We applied the resummation formalism to the fully- 
differential cross section for the production of the Higgs 
boson at the LHC, combining the most advanced per- 
turbative information available at present: soft-gluon 
resummation at NNLL accuracy and fixed-order pre- 
diction at NLO QCD. The numerical results show a 
sizeable resummation effect at intermediate qx, a qx- 
shape mildly dependent on the rapidity of the Higgs 
boson and an overall stability of the cross section with 
respect to scale variation and to the inclusion of higher 
perturbative orders. 

Acknovifledgements. This work was supported 
by the Deutsche Forschungsgemeinschaft under SFB 
TR-9 "Computergestiitzte Theoretische Teilchenphysik" 



References 

1. For a review on Higgs physics in and beyond the 
Standard Model, see J. F. Cunion, H. E. Haber, 
G. L. Kane and S. Dawson, The Higgs Hunter's 
Guide (Addison- Wesley, Reading, Mass., 1990); 
M. Carena and H. E. Haber, Prog. Part. Nucl. 
Phys. 50 (2003) 63; A. Djouadi, report LPT-ORSAY- 
05-17 [ hep-ph/0503172 , report LPT-ORSAY-05-18 
hep-ph/0503173i . 



Colliders - Higgs Phenomenology 



Contributed Talk 



2. ATLAS Coll., ATLAS Detector and Physics Per- 
formance: Technical Design Report, Vol. 2, report 
CERN/LHCC/99-15 (1999); S. Asai et al., Eur. Phys. 
J. C 32S2 (2004) 19; CMS Coll., CMS Physics Tech- 
nical Design Report: Physics Performance, Vol. 2, re- 
port CERN/LHCC/2006-021 (2006). 

3. M. Carena et al.. Report of the Tevatron Higgs work- 
ing group, liep-pli/0010338; CDF and DO Collabora- 
tions, Results of the Tevatron Higgs Sensitivity Study, 
report FERMILAB-PUB-03/320-E; V. M. Abazov 
et al. [DO Coll.], Phys. Rev. Lett. 96 (2006) 011801; 
A. Abulencia et al. [CDF Coll.], Phys. Rev. Lett. 97 
(2006) 081802; The TEVNPH working group [for the 
CDF and DO Collaborations], hep-ex/0612044 

4. R. K. Ellis, I. Hinchliffe, M. Soldate and J. J. van 
der Bij, Nucl. Phys. B 297 (1988) 221; U. Baur and 
E. W. Glover, Nucl. Phys. B 339 (1990) 38. 

5. D. de Florian, M. Grazzini and Z. Kunszt, Phys. Rev. 
Lett. 82 (1999) 5209. 

6. V. Ravindran, J. Smith and W. L. Van Neerven, Nucl. 
Phys. B 634 (2002) 247. 

7. C. J. Glosser and C. R. Schmidt, JHEP 0212 (2002) 
016. 

8. V. Del Duca, W. Kilgore, C. Oleari, C. Schmidt and 
D. Zeppenfeld, Nucl. Phys. B 616 (2001) 367, Phys. 
Rev. D 67 (2003) 073003. 

9. J. Smith and W. L. van Neerven, Nucl. Phys. B 720 
(2005) 182. 

10. C. Anastasiou, K. Melnikov and F. Petriello, Phys. 
Rev. Lett. 93 (2004) 262002, Nucl. Phys. B 724 
(2005) 197. 

11. S. Catani and M. Grazzini, Phys. Rev. Lett. 98 (2007) 
222002. 

12. Y. L. Dokshitzer, D. Diakonov and S. I. Troian, Phys. 
Lett. B 79 (1978) 269, Phys. Rep. 58 (1980) 269. 

13. G. Parisi and R. Petronzio, Nucl. Phys. B 154 (1979) 
427. 

14. G. Curci, M. Greco and Y. Srivastava, Nucl. Phys. B 
159 (1979) 451. 

15. J. C. Collins and D. E. Soper, Nucl. Phys. B 193 

(1981) 381 [Erratum-ibid. B 213 (1983) 545], Nucl. 
Phys. B 197 (1982) 446. 

16. J. Kodaira and L. Trentadue, Phys. Lett. B 112 

(1982) 66, report SLAC-PUB-2934 (1982), Phys. 
Lett. B 123 (1983) 335. 

17. J. C. Collins, D. E. Soper and G. Sterman, Nucl. 
Phys. B 250 (1985) 199. 

18. S. Catani, E. D'Emilio and L. Trentadue, Phys. Lett. 
B 211 (1988) 335. 

19. S. Catani, D. de Florian and M. Grazzini, Nucl. Phys. 
B 596 (2001) 299. 

20. D. de Florian and M. Grazzini, Phys. Rev. Lett. 85 
(2000) 4678, Nucl. Phys. B 616 (2001) 247. 

21. G. Bozzi, S. Catani, D. de Florian and M. Grazzini, 
Phys. Lett. B 564 (2003) 65. 

22. G. Bozzi, S. Catani, D. de Florian and M. Grazzini, 
Nucl. Phys. B 737 (2006) J73^ 

23. http:/ /theory.fi. infn.it/grazzini/codes. html" 

24. G. Bozzi, S. Catani, D. de Florian and M. Grazzini, 
arXiv:0705.3887 [hep-ph]. 

25. A. D. Martin, R. G. Roberts, W. J. Stirling and 
R. S. Thorne, Phys. Lett. B 604 (2004) 61. 



